ENHANCED DELAY TO BIFURCATION 
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Abstract. We present an example of slow- fast system which displays a full 
open set of initial data so that the corresponding orbit has the property that 
given any e and T, it remains to a distance less than e from a repulsive part 
of the fast dynamics and for a time larger than T. This example shows that 
the common representation of generic fast-slow systems where general orbits 
are pieces of slow motions near the attractive parts of the critical manifold 
intertwined by fast motions is false. Such a description is indeed based on the 
condition that the singularities of the critical set are folds. In our example, 
these singularities are transcritical. 



A first approximation for the time evolution of fast-slow systems dynamics is 
often seen as follows. A generic orbit jumps close to an attractive part of the 
equilibrium set of the fast dynamics. It evolves then slowly close to this attrac- 
tive part until that, under the influence of the slow dynamics, this attractive 
part bifurcates into a repulsive one. Then the generic orbit jumps to another at- 
tractive part of the fast equilibrium set until it also looses its stability and then 
jumps again in search of another attractor and so on. This approach is indeed 
only a first approximation. It is a quite meaningful one because in this setting 
one can explain many phenomena like for instance hysteresis cycles, relaxation 
oscillations, bursting oscillations and more complicated alternance of pulsatile 
and surge patterns of coupled neurons populations. It can also serve as a good 
representation model of the so-called ultrastability introduced in Cybernetics by 
R. Ashby ([]]) (as the property that some systems display an adaptation under 
some external influences to switch in a reversible or non-reversible manner from 
a stable state to another one) and used in life sciences by H. Atlan (j2], [10]). 

Also, this approach does not take into account canards. Discovered by E. 
Benoit, Callot, M. and F. Diener ([3]) using non-standard analysis, canards are 
limit periodic sets (limit when e — > 0) of the van der Pol system: 
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(1.2) y = X -c(e) 

where c(e) ranges between some bounds: 

(1.3) c + exp(-a/e) < c(e) < c + exp(— /3/e), (a > (3 > 0). 

The surprising aspect is that a part of the canard coincides with a repulsive piece 
of the critical manifold. More recent contributions of Dumortier and Roussarie 
(0) yield new (standard) proofs for the existence of such limit-periodic sets. 
There are now several evidences showing the relevance of this notion to explain 
experimental facts observed in physiology. Note that canards occur in two- 
parameters families in some narrow set of parameters. 

The system we present here, although quite simple, displays a kind of recurrence 
of the canard effect which enhances the delay to bifurcating. 

2. Transcritical Dynamical Bifurcation 

The classical transcritical bifurcation occurs when the parameter A in the 
equation: 



(2.4) x = -Xx + x 2 , 

crosses A = 0. Equation 12.41 displays two equilibria, x = and x = A. For 
A > 0, x = is stable and x = A is unstable. After the bifurcation, A < 0, x = 
is stable and x = A is unstable. The two axis have "exchanged" their stability. 
The terminology "Dynamical Bifurcation" refers to the situation where the bi- 
furcation parameter is replaced by a slowly varying variable. In the case of the 
transcritical bifurcation, this yields: 



(2.5) x = —yx + x 2 

(2.6) y = -e, 

where e is assumed to be small. 
This yields 

x = -(-et + y )x + x 2 , (y = y(0)) 
which is an integrable equation of Bernoulli type. Its solution displays: 

x(t) = — ^M~nt)j {xo = x{0) 

1 — Xq / exp[— Y (u)\du 

ft r 1 t 2 

Y(t) = / y(s)ds = / (-es + y )ds = -e— + y t. 
Jo Jo 2 

If we fix an initial data (x ,y ), y > 0, < x < yo/2, and we consider the 

solution with this initial data we find easily that it takes time t = y^jt to reach 

the axis y = 0. If xq is quite small, that means the orbit stays closer and closer 

of the attractive part of the critical manifold untill it reaches the axis x = y and 

then coordinate x start increasing. But now consider time cyo/e, 1 > c > 2. 
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Then a straightforward computation shows that Y(t) = c(l — §) — = -, and 
that 

f2yo/e t 2 , , 2 

o / exp(e— - y *)<« = — ■ 
Je 2 y 

This shows that if e is small enough and if the initial data satifies Xq < yo/2, 
then: 

g exp(-fc/e) 
1 - 2zo/yo 

This yields that, despite the repulsiveness of the axis x = 0,y < 0, the orbit 
remains for a very long time close to x = 0, indeed < xo- Note that after 
a larger time the orbit blows away from this repulsive axis. This phenomenon, 
although quite simply explained, is of the same nature as the delay to bifurcation 
discovered for the dynamical Hopf bifurcation. See for instance ([3], [6], [9], 
[llj). This well-known effect is instrumental in the example we construct in 
this article. Some related work has been done in computing exit points through 
the passage of single turning points, for example [I], [TJ. 



3. AN EXAMPLE OF SYSTEM WITH ENHANCED DELAY 
Consider the equation: 

x — (1 — x 2 )(x — y) 
y = ex. 

The fast dynamics displays the invariant lines x — — 1, x — 1, and y = x. 
A quick analysis shows that, as the slow variable y varies, the fast system 
undergoes two transcritical bifurcations near the points (—1,-1) and (1,1). As 
we recalled in the first paragraph, a typical orbit near (—1, —1) first displays a 
"delay" along the repulsive part (x = —l,y < —1) of the slow manifold. Then, 
by hysteresis, it jumps to the attractive part (x — l,y < 1) till it reaches the 
other transcritical bifurcation where it again displays another delay along the 
repulsive part (x — 1, y > 1). Then it jumps again to (x = —1, y > 1) and starts 
again. There is such a mechanism of successive enhancements of the delay after 
several turns generated by the hysteresis. After this intuitive explanation, we 
give now a formal proof of the: 
Theorem 

For all initial data inside the strip — 1 < x < 1, for all S and for all T, the 
corresponding orbit spends a time larger than T within a distance less than S 
to the repulsive part of the slow manifold. 
Proof 

Inside the strip | x \< 1, it is convenient to use the variable u: x = tanhu. The 
system yields the equations: 

u = tanhu — y 



y = etanhu. 
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Consider the function 

$(u, y) = -(tanhu — y) 2 + em(coshu), 
and its time derivation along the flow. This displays: 



dt ^ ' ^ ^ coshu ^ 



Hence the function $ is strictly increasing along the flow (Lyapunov function 
for the flow). 

Note now that if (u(t),y(t)) is a solution then (— u(t), —y(t)) is also a solution. 
To study the orbits of the system, we can restrict to initial data u = uo > 
and y = y . The first step of the proof is to show that all orbits intersect both 
axes u = and y — u in infinitely many points. 
Assume first y > tanhu . Then -u(O) < 0. But 



shows that y— tanhu grows hence remains positive. Assume that u would remain 
always positive. Then, as y — tanhu > 0, u is monotone decreasing. Hence there 
exists I such that u — > I as t — > +oo. 

If I > 0, y — etanhw implies (via the mean value theorem) y — > +oo but then 
ii — > — oo and (mean value theorem) contradiction with u > 0. 
If I — 0, u + y — > 0. But y is monotone increasing. If y — > +oo, then -u — > — oo 
and again contradiction. If y tends to a finite limit m, then ii — > — m and again 
contradiction. Hence all orbits with initial data (uo,yo) with y > tanh-u > 
intersect the axe u — 0. 

Consider now the case y < tanh-u - The variable u is first strictly increasing (as 
soon as y < tanhu), hence positive and so y is increasing. Assume that tanhu— y 
would remain positive along the orbit. Then as t — > +oo, u would tend to a 
limit m (eventually m = +oo).As y — > m, mean value theorem would imply 
y — > +oo and again a contradiction with -u — > — oo. So the orbit necessarily 
intersects the axe y = tanhu and ultimately the axe u = by the preceding 
argument. By symmetry, we can also show the existence of two sequences of 
times (t n ) and (9 n ) such that: 



In the second part of the proof we show that the sequence (a n ) is unbounded. 
Consider now the function 

2$(-u, y) = w(x, y) = (x — y) 2 — eln | 1 — x 2 |, 

which satisfies: 



d_ 

dt 



(y — tanhu) = etanh-u + 



y — tanhu 




w = 2(1 - x 2 ){x - yf / 1. 



ENHANCED DELAY TO BIFURCATION 



5 



As w is strictly increasing, this yields: 

w(t n ) =a 2 n < w(9 n ) = -eln(l - < w(t n+1 ) = a 2 n+1 , 
with £ n = x(0 n ). Integration along the flow of (w) = 2x(x — y) yields: 

ftn+l rtn+1 

= / 2x(x - y)dt = 2 xydt 

ftn+l 

= 2e x 2 dt < 2e(t n -t n+1 ). 
j t n 

This shows that if (t n ) converges to a finite limit then so does (a„) and (£ n ). 
Now integration along the flow of y = eu yields: 

rtn+1 

a n + a n +i < e / I tanhw | dt < e(t n+1 - t n ), 

and this shows that the sequence of times (t n ) is necessarily unbounded. Now 
assume that the sequence (a„) would be bounded. Then, as the sequence (t n ) 
tends to +00, the function w would be bounded on the orbit. But then so would 
be both the two functions (x — y) 2 and — ln(l — x 2 ). But then there would exist 
a constant a such that (1 — x 2 ) > a along the orbit and 

w > 2a(x — y) 2 > 2aw, 

hence e~ 2at w(t) increasing and contradiction with the fact that w would be 
bounded. 

Last step is classical in slow-fast dynamics. Consider the unbounded sequence 
of points (0, a n ) on the orbit. By Tikhonov's theorem there is an invariant curve 
which passes near that point. By Takens's theorem the flow is conjugated to 
the fast flow along this invariant curve. Hence necessarily there is for all orbits 
inside the strip, all 5 and all T a part of the orbit which remains at a distance 
less than 5 of the repulsive parts of the boundary of the strip for a time larger 
than T. 

4. ASYMPTOTICS OF THE SYSTEM OUTSIDE THE STRIP 

Theorem 

Given any initial point (x ,y ) outside the strip | x \< 1, the corresponding 
orbit is asymptotic to y — x. 
Proof 

We can always assume that x > 1 because the system is symmetric relatively 
to the origin. The equations yield: 

e^ = (x- -)(y-x). 
dy x 

So if y > x and (xo, yo), (xi,yi) are two points on the same orbit with xq < x±, 
we get: 

1 dx 1 

(ar - — ){y ~ x) < e— < ( Xl - — ){y - x). 
x ay Xi 
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Set: 



this yields 



hence 



ttj = -(Xi - — ), 2 = 0,1 
€ Xi 

dx 

— - a (y - x) > 0, 
^ - ai (y-x) < 0, 

^-(e aoy x) > a ye a ° y 
dy 

-^-(e aiy x) < ai ye aiy . 
dy 

Integration between y and y\ yields 

e aoyo (yo -x -—)> e a ^( Vl - Xl - — ), 

and 

e aiyo (Vo -x -—)< e a ^( Vl - Xl - — ). 

The second inequality shows that if y = x , then 

e aiyo >e Qm (l-o;i(yi-a;i)), 

and thus that the orbit stays above the line y = x. If we now choose x Q (and 
ao) large enough, the first inequality displays: 

1 

[yi -xi) < —■ 
2a 

This shows that the orbit is asymptotic to y — x. 

References 

[1] R. Ashby, Cybernetics, New York, Academic Press, 1959. 

[2] H. Atlan, L 'organisation biologique et la theorie de I 'information, Paris, Hermann, 1972. 
[3] S.M. Baer, T. Erneux and J. Rinzcl, The slow passage through a Hopf bifurcation: Delay 

memory effects and resonance SIAM Journal on Applied mathematics, vol. 49, nl (Feb 

1989), 55-71. 

[4] E. Benoit, Equations differentielles: relation entree-sortie Comptes rendus de l'Acadcmic 

des Sciences dc Paris, 293, serie I, 293-296 (1981). 
[5] E. Benoit, J.-L. Callot, F. Diener and M. Diener, Chasse au canard. Collectanea Math- 

ematica, Barcelone, 31(1-3); 37-119 (1981). 
[6] B. Candelpcrghcr, F. Diener and M. Diener, Retard a la bifurcation : du local au global hi 

Bifurcations of Planar Vector Fields proceedings of the Luminy conference, JP Frangoise 

and R. Roussarie (Edts), Lecture Notes in Mathematics 1455, Springer, Berlin (1990), 

1-19. 

[7] P. De Maesschalk, F. Dumortier, Time analysis and entry-exit relation near planar turn- 
ing points. Journal of Diff. Equations, vol. 215, Issue 2, 225-267 (2005). 

[8] F. Dumortier and R. Roussarie, Canard cycles and center manifolds with an appendix 
by Cheng Zhi Li. Mem. Amer. Math. Soc. 121 (1996), n577. 



ENHANCED DELAY TO BIFURCATION 



7 



[9] T. Erneux, E.L. Reiss, L.J. Holden, M. Georgiou, Slow passage through bifurcation and 
limit points. Asymptotic theory and applications. Dynamic Bifurcations (Luminy, 1990), 
Ed. E. Benoit, Lecture Notes in Mathematics 1493, Springer, Berlin (1991), 14-28. 

[10] J. -P. Frangoise, Oscillations en biologic : Analyse qualitative et modeles, Collection: 
Mathematiques ct Applications 46, Springer, (2005). 

[11] H. Zoladek, Remarks on the delay of the loss of stability of systems with changing param- 
eters in Bifurcations of Planar Vector Fields proceedings of the Luminy conference, JP 
Frangoise and R. Roussarie (Edts), Lecture Notes in Mathematics 1455, (1990), 393-396. 

Laboratoire J.-L. Lions, UMR 7598, CNRS, Universite P.-M. Curie, Paris6, 
Paris, France 

E-mail address: Jean-Pierre.Francoise@upmc.fr 



